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I . Abstract 
^ ■ 

\ In this paper we study the hnear stabihty of the relative equiUbria for 

^ ' homogeneous and quasihomogeneous potentials. Firstly, in the case 
'— the potential is a homogeneous function of degree —a, we find that any 

\ relative equilibrium of the n-body problem with a > 2 is spectrally 

J> ' unstable. We also find a similar condition in the quasihomogeneous 

^ ■ case. Then we consider the case of three bodies and we study the 

Q>^ , stability of the equilateral triangle relative equilibria. In the case of 

tJ" ' homogeneous potentials we recover the classical result obtained by 

Q> ■ Routh in a simpler way. In the case of quasihomogeneous potentials 

, we find a generalization of Routh inequality and we show that, for 

f-*) \ certain values of the masses, the stability of the relative equilibria 

' depends on the size of the configuration. 
• ^ . 

^ \ Keywords: Three body problem, central configurations, linear stability, 

■ homogeneous potentials, quasihomogeneous potentials. 



1 Introduction 

The study of the stability of the relative equilibria is very important in the 
analysis of Hamiltonian systems. In the last few decades some general meth- 
ods have been developed to tackle this problem, namely the energy- Casimir 
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method and the energy momentum method. The energy-Casimir method is 
an adaptation of the classical Lagrange-Dirichlet method [TT] that was devel- 
oped by Arnold fl\ and applied to analyze the stability of stationary flows of 
perfect fluids. Such method was also applied to a variety of other problems. 
However this technique relies heavily on a supply of Casimir functions. The 
energy-momentum method was developed to overcome this difficulty, it is 
very powerful and was applied successfully to many problems. 

Unfortunately there are several problems that cannot be tackled using the 
energy momentum method. Indeed in some cases the second variation of the 
energy is indeterminate. In particular it is known that in celestial mechanics 
the stability of the relative equilibria cannot be studied using the energy 
momentum method. However one can study the linear or spectral stability. 
The actual stability of relative equilibria is a major open problem in celestial 
mechanics. Several authors (see for example |12j) think that all the relative 
equilibria are nonlinearly unstable because Arnold's diffusion is believed to 
be a feature of the n-hodj problem. Remarkably it is very difficult to study 
Arnold's diffusion numerically and the existence of Arnold's diffusion in the 
n-body problem has not been proved beyond doubt. 

In this paper we consider the problem of n bodies interacting by means 
of homogeneous and quasihomogeneous potentials. In the case the potential 
is a homogeneous function of degree —a we find that any relative equilibrium 
is unstable if a > 2. In the case of quasihomogeneous potentials we find a 
similar, but more complex condition. Then we study the equilateral triangle 
relative equilibria and their stability. The energy-momentum method fails 
to provide informations concerning the stability of the relative equilibria 
since the second variation of the amended potential is indefinite. Therefore 
we study the linear stability of those solutions. The linear stability of the 
Lagrangian triangle solution of the Newtonian three body problem was first 
studied by Gascheau in 1843 Later Routh |14| studied the linear stability 
of the same solutions in the case of homogeneous potentials. He proved that 
linear stability is achieved only when a < 2 and the masses satisfy 

mim2 + minis + 17121123 ^ 1 /2 - a V 

where the potential is a homogeneous function of degree —a with a > 0. In 
this work, using a considerably simpler derivation based on an idea of Moeckel 
(see [H]), we recover the classical result of Gascheau and Routh (for spec- 
tral stability), and we analyze it. In particular we find that the equilateral 
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triangle relative equilibria can be spectrally stable for < a < 2 if one of 
the masses is much larger than the other two. An interesting feature appears 
when < a < 14 — 8\/3, indeed, in this interval, there is spectral stability 
not only in presence of a dominant mass but also when one of the masses is 
sufficiently small. This is a somewhat surprising result that stresses the dif- 
ference between the Newtonian potential and other homogeneous potentials. 
Indeed it is widely believed that, in the Newtonian Moeckel conjec- 

tured yjj, a relative equilibrium can be linearly stable only if it contains a 
mass significantly larger than the other masses. 

Then we use similar methods to generalize Routh condition to the case 
of quasihomogeneous potentials. Quasihomogeneous potentials ([SI [ID]) have 
been widely studied in recent years since they are suitable to describe a 
number of phenomena. Among the other we recall some of the most rep- 
resentative: the so called Manev (see [H El [3], [9] and references therein) 
and Schwarzschild (see [HI [17]) potentials. The Manev potential and the 
Schwarzschild potentials are particularly interesting because, under suitable 
conditions, they describe the motion of a test particle in a gravitational field. 
In particular it was shown that the Manev potential can describe the preces- 
sion of the perihelion of Mercury with the same accuracy as general relativity 
(see [8]). The Schwarzschild potential can also be used for similar calculations 
(for example see [T7]). 

The Lennard- Jones potential ([2]) is also a very interesting potential used 
in molecular dynamics to simulate many particle systems as for example 
solids, liquids and gases. The relative equilibria for this potential were studied 
in [2] and will not be considered here. 

In the case of quasihomogeneous potentials we find that spectral stability 
is achieved only when 



(a2 - 2ay+^ + {}? - 26)r^+2 



< 



and ^ 

(mi + m2 + m3)2 " 3 \h{h + 2)r^^-^ + a{a + 2) ) ^' 
where Tq is the mutual distance between any two bodies and the potential 
is the sum of a homogeneous function of degree —a and of one of degree —h 
with a > 6 > 0. We denote with / the left hand side of inequality ([2]). 

Analyzing the conditions above we find that, for certain values of the 
masses, the stability of the Lagrangian triangle solutions depend on the size 
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of the triangle. In particular we show that if a > 6 > 2 the triangle solutions 
are unstable for any value of the masses. Ifa>2>6>0 the above solutions 
are unstable when / > l/3[(6-2)/(6+2)]2 and when / < l/3[(6-2)/(6+2)]2 
and tq < zl (for some zl G M"*"). Moreover they are spectrally stable when 
/ < l/3[(6-2)/(6+2)]2 and ro> z^. If 2 > a > 6 > there are three different 
cases depending on the value of the masses. When / > l/3[(6 — 2)/(6 + 2)]^ 
the solutions are unstable, while when / < l/3[(a — 2)/(a + 2)]^ they are 
spectrally stable. When l/3[(a-2)/(a + 2)]2 < / < l/3[(6-2)/(6 + 2)]2 they 
are unstable for ro < tq*, but spectrally stable for tq > Tq for some Tq G M+. 

The study of the stability of relative equilibria for quasihomogeneous 
potentials is quite interesting. Indeed the Trojan asteroids are found at the 
Lagrangian point of the Sun- Jupiter system. The gravitational interaction is 
usually considered to be the Newtonian one in most of the studies concerning 
the Trojans. However if one wants to consider the problem in the framework 
of General Relativity (see [15] for a study of the three-body problem in the 
Post-Newtonian approximation of General Relativity ) then one might get 
some surprising results. Since the Schwarzschild and the Manev potentials 
are used to approximate general relativity the results presented in this paper 
suggest that it might be impossible to find, in nature, equilateral triangle 
relative equilibria smaller than a certain size. This is because the relative 
equihbria studied in this work (in the case of quasihomogeneous potentials) 
are unstable if the interacting bodies are too close to each others. 

This paper is organized as follows. In the next section we write the 
equation of motion in Cartesian coordinates and in uniformly rotating co- 
ordinates. We define the relative equilibria for the n-body problem and we 
formulate the conditions for their spectral stability. In section 3 we study 
the case of homogeneous potentials. First we consider the n-hodj problem 
in general and then for n = 3 we recover the inequality of Routh to describe 
the stability of the Lagrange triangle relative equilibria. In the last section 
we consider quasihomogeneous potentials. We first study the n-body prob- 
lem and then we find a condition for the spectral stability of the Lagrange 
triangle solutions and we discuss it. 
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2 Relative Equilibria and Their Stability 



2.1 Relative Equilibria 

We let the mass and position of the n bodies be given by rrii G M"*" and 
qj G M^, where z = 1, . . . , n. Let q = (qi, . . . , q„) G M^". Newton's equations 
for the zth body are 

dU 

mifii = — 

where V({q) is the potential function. We let the momentum of each body be 
Pi = rriiqi and let p = (pi, . . . , p„) G M^". The equations of motion can be 
written as 



q = M~^p -- 
p = VU{q) 



OH 

dp 



OH 
dq 



(3) 



where M is the 2n x 2n mass matrix diag(mi, mi, . . . , m„, m„), V denotes 
the Euclidean gradient in M^*^ and H{q,p) is the Hamiltonian function 

n II ||2 



1=1 



Let J and R{0) denote the 2n x 2n block diagonal matrices with n equal 
2x2 blocks of the form 



K 



1 
-1 



and e^' 



cos sin 
— sin 9 cos 9 



respectively. To introduce coordinates that uniformly rotate with constant 
angular velocity a), we let Xj = (\ie^^^ and = pjc'^^'*. This is a symplectic 
change of variable and thus preserves the Hamiltonian structure of system 
([3]). The equations of motion in the new variables are 



X = CjJx + M = 
y = VU{x) + i2}Jy 



dH 

dy 
dH 
dx 



(4) 



5 



where x(t) = R{(ljt)q{t), y(t) = R{Cjt)p{t) and H{x,y) is the Hamiltonian 
function: 

H{x, y) = -y^M-'y - U{x) - tox^Jy. 

An equihbrium point of system (jlj) corresponds to a periodic solution in the 
n-body problem consisting of a configuration of masses which rotates rigidly 
about its center of mass. An equilibrium (x, y) of system (jlj) must satisfy 
y = —CjMJx and 

VU{x) + u'^Mx = (5) 

A relative equilibrium is a configuration x G M^" which satisfies the algebraic 
equation ([5]) for some value of u. This can be viewed as the condition for a 
critical point of the restriction of the potential to {x G M^" : x'^Mx = c} for 
any constant c, where u'^ plays the role of a Lagrange multiplier. 

In the following we will consider the potential U to be either a homoge- 
neous potential of the form 

and 

where < 6 < a, or a quasihomogeneous potential defined as a sum of the 
homogeneous potentials above: 

U{q) = Viq) + W{q). (8) 

A configuration a; is a relative equilibrium for the quasihomogeneous potential 
U provided that 

\/U{x) = -uj^Mx (9) 
for some constant uj, while it is a relative equilibrium for V and W if 

\/V{x) = -ujIMx and \/W{x) = -ujIMx (10) 

respectively. Since the potentials V and W are homogeneous of degree —a 
and —h respectively the constants uji and UJ2 are determined by: 

2_aV{x) 2_,W{x) 
- x^ - ^x^ 
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Furthermore if x is simultaneously a relative equilibrium for V, W and U 
then 

VU{x) = VV{x) + VW{x) = -{ujI + uol)Mx 

and 

2 2,2 
U = UJi+ UJ2- 

2.2 Linear and Spectral Stability 

A relative equilibrium x is linearly stable if the origin is a stable solution of 
the linearization at x of system (jlj). A necessary condition for x to be linearly 
stable is that all the eigenvalues of the linearization 



S 



DWU CjJ 



(12) 



are either zero or purely imaginary. This weaker condition is called spectral 
stability. 

The characteristic polynomial of S, -P(A), is of degree An and is an even 
polynomial, since S" is a Hamiltonian matrix. Let v be an eigenvector of S 
with eigenvalue A and write v = {vi, V2), where vi,V2 G C^". The eigenvector 
equation reduces to 

V2 = M{XI - LUJ)VI 

Avi = ^ ' 

where 

A= M-^DVU + {iu'^ - X"^)! + 2XiuJ (14) 

and I is the identity matrix. Consequently, to obtain the eigenvalues of 
S, one needs only take the determinant of A and find the roots, namely 
P(A) = det(^). Following Moeckel we introduce the normalized eigenvalues 
fi = X/\uj\. These satisfy the equation 

det(M-^DVW + (l-/i2)/ + 2/iJ) = 0. (15) 

Note that the stability condition is unchanged since the normalization factor 
is real and positive. 

Two vectors v and w are called M-orthogonal if v'^MW = 0. One 
can show that J is antisymmetric and M^^DVU is symmetric in any M- 
orthonormal basis. Using an M-orthogonal basis and taking the transpose 
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of A does not change the determinant (fT3]) and thus the determinant is an 
even function of fi. Let z = ^"^ and let G{z) be the polynomial of degree 2n 

Gifi^) = det{M-^DVU + (1 - /i^)/ + 2^J) = 0. (16) 

Then x is spectrally stable if and only if all the roots of G{z) are either zero 
or real and negative. G{z) is called the stability polynomial. 

In order to study spectral stability it is convenient to obtain factorizations 
of the stability polynomial G{z). Moeckel's idea is that such a factorization 
can be obtained by finding subspaces of M^" which are simultaneously invari- 
ant for J and M'^DVU: 

Proposition 1. Suppose that W G M^" is a invariant suhspace for both J 
and M~~^DVl/( . Then the stability polynomial can be factored into two even 
polynomials in G{z) = Gi{z)G2{z) where Gi{z) and G2{z) are given by l[TS\) 
with the matrices involved restricted to the subspaces W and W-^ = {f G 
R2" : v^Mw = 0y w eW}. 

The invariant subspaces W must have even dimension. The simplest case 
is dimension two. In this paper we will only need to use invariant subspaces 
of dimension two. In this case we have the following 

Proposition 2. Suppose W G M^" is a two-dimensional subspace which is 
simultaneously invariant for J and M^^DVU . Let the eigenvalues of the 
restriction of lj~'^M~^D'VU be rj and ^. Then G{z) has a quadratic factor 

Q(z) = z^ + az + l3 

where a = 2 — rj — ^ and (5 = {1 + ri){l + ^) . The corresponding roots are real 
and negative if and only if 

a>0 (3>0 - 4/? > 0. 



3 Homogeneous Potentials 
3.1 General Case 

Consider the potential V{x). Then the matrix DVV"(x) is of the form: 



A 



An 



Air. 
A 
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where Aij is the 2x2 matrix given by 



^kk — ~ / ^ Ajk 




[/ - (a + 2)ujkuji.] if j ^ k 



(17) 



where u^^ = ^^^-^ and djk = ||xj — Xfc||. 

Using the fact that the diagonal blocks of A are the negative of the sum of 
the blocks in the corresponding rows it is clear that both v = (1, 0, . . . , 1, 0) 
and w = (0, 1, . . . , 0, 1) are in the kernel of A. Therefore spanjf , w} is a 
two-dimensional invariant subspace for both M^^A and J. The eigenvalues 
of the restriction of ujfM'^A are rj = ^ = 0. Applying Proposition [2] one finds 
that the the roots of Q{z) are real and negative. Substituting rj = = into 
Q{z) yields the repeated root z = —1. The corresponding eigenvalues are 
fi = ±i. These values are a result of a drift in the center of mass, and are 
present in any relative equilibrium. 

Another two-dimensional invariant subspace comes from the configuration 
itself. Consider v = x and w = Jx. Equation (ITU]) with the homogeneity of 
the potential gives: 



so V is an eigenvector of M~^B with eigenvalue (a + 1)^1. On the other 
hand w = (i?~^(0))'x where R{0) is the 2n-dimensional rotation operator 
introduced above. Moreover 



because of the rotation invariance of the potential. Differentiating at 6 = 



and thus w is an eigenvector with eigenvalue —uf. Choosing rj = (a + 1) and 
= — 1 we find that a = 2 — a, and the roots of Q{z) are positive if a > 2. 
This proves the following 

Theorem 1. Any relative equilibrium ofV{x) is spectrally unstable if a> 2. 



M-^DW{x)v 



-{a + l)M'^VV{x) 

(a + 1)M~^ujIMx = {a + 1)ujIv, 



(18) 



VV{R-^{e)x) = R-\9)W{x) 



gives: 



M-^DVV{x)w = M-^JVV{x) = -M-^JuolMx = -ufw 



(19) 
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Figure 1: (a) Some level curves /(mi, m2, m^) = c in the triangle of normal- 
ized masses. Regions of stability (shaded) and instabihty in the triangle of 
normalized masses for the homogeneous potential V with 14 — 8\/3 < a < 2. 
(b) Regions of stability (shaded) and instability for the potential V with 
< a < 14 - 8\/3. 



3.2 Lagrangian Triangle Relative Equilibria 

Now let n = 3. Then the equilateral triangle configurations with the center 
of mass at the origin are relative equilibria, for any choice of the masses 
mi, 7712 and 7723. Moreover, since the potential is homogeneous, the stability 
analysis is not affected by the size of the configuration. Therefore, without 
loss of generality, we can fix the mutual distances. Let djk = Vs. Then 
the potential energy is V{x) = (mim2 + minis + m2m3) / {y/S)"" , the moment 
of inertia is x^Mx — 3(mim2 + mim^ + 777.2/77.3)//^, and uf — a/i/(-\/3)"'^^. 
Assume that the particles are at positions Xj — Xj — c where c is the center 
of mass, and 

x[^{l,0) 4=(-l/2,v^/2) 4 = (-1/2, -v^/2) 

Since the formulas for A are translation invariant the x'jj, can be used instead 
of the Xjk and one finds that C = ijJi'^M~^A is: 





C12 


Ci3 


C21 


C22 


C23 


C31 


C32 


^33 
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where 



Cii = 



Cl3 = 



C22 - 



C31 = 



4m 
1 

1 

4m 
1 

4m 



(2 + 3a)(m2 + ms) v/3(a + 2) (m3 - m2) 
v'3(a + 2)(m3 — 1712) —(2 — o){m2 + ma) 
-(2 + 3a)m3 -^{a + 2)7713 

— \/3{a + 2)m3 (2 — 0)7773 
(2 + 3a)777i — 47773 -v^{a + 2) mi 

— y/3{a + 2)mi (a — 2)mi + 4(1 + a)m3 

-(2 + 3a)mi -V^(a4-2) 7771 

— \/l3(a + 2)mi (2 — a)777i 



(^12 — 



4m 



C21 — 



4m 



^"23 ■ 



4m 



C32 — 



4m 



-(2 + 3a)7n2 ^3(0 + 2)7712 
v^(a + 2)7772 (2 — c^)7772 

-{2 + 3a)777i y3(a + 2)7771 

V^(ci + 2)7711 (2 — a)mi 

47713 
-4(1 + 0)7713 

47772 

-4(1 + a)77i2 



and 



C33 — ■ 



4m 



(2 + 3a)777i - 47772 V3(a + 2)mi 

v/3(a + 2)7771 (a - 2)7771 + 4(1 + a)7772 



Four eigenvalues of C are 0, 0, —1, a + 1 that we found in subsection 3.1. 
Let us denote with rj and ^ the remaining eigenvalues. The eigenvalues 7] and 
^ can be determined, however we will only need to know the value of ?7 + ^ 
and T]^. The trace ofCisa + 1 — l+?7 + ^ = 2a, i.e. r] + ^ = a, thus 



a 



2-^-r] = 2-a. 



To find /5 = (?7 + 1)(^ + 1) = r]^ + 1 + a we need to compute rj^. 

The sum of all possible product of pairs of eigenvalues is — a — 1 + rj^. 
This sum can also be found as one-half the difference of the square of the 
trace of the matrix above and the trace of its square. This computation 
yields 



and hence 



[(lla^ — 4a — 4)(mim2 + 77117713 + 7722^3) + (4a^ — 8a — 8) {ml + m\ + 772,3)], 

?) {a^ + Aa + A){mim2 + mim^ + m2rn^) 
^"4 'i^^ ■ 

Since a > and /? > the equilateral triangle is spectrally stable if and only 
if - 4/3 > 0, that is 

777 1 7772 + 777x7773 + 7772^3 ^ 1 / 2 - a \ ^ 

(7771 + 7772 + 7773)2 -3V2 + ay ■ 

Let f {mi, 1712,7713) be the left hand side of equation fl20l) . Then the level 
curves f{mi, 7772, 7/7-3) = c can be plotted on the triangle of normalized masses. 
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see Figure 1(a) The level curves for levels less than 1/4 consist of three 



curves (arcs of circle) near the corners of the triangle, while the levels for 
c > 1/4 determine a single closed curve, a circle centered in the orthocenter 
of the equilateral triangle. The value 1/4 corresponds to a = 14 — 8v^ ^ 
0.14359354. 



Figure 1(a), that depicts the triangle of normalized masses for the ho- 
mogeneous potential V when 14 — 8\/3 < a < 2, can be interpreted in the 
following way . The shaded region represents the region of stability, i.e. the 
values of the masses for which the triangle solutions are stable. Thus, when 
14-8^3 < a < 2, the Lagrangian triangle solutions are stable only for the 
mass triplets in the corners of the triangle, i.e. when one mass is considerably 
larger than the other two. In the Newtonian case, namely when a = 1, they 
are unstable for the majority of mass triplets. 



Figure 1(b) depicts the triangle of normalized masses when < a < 



14 — 8\/3. As before the shaded region represents the region of stability. 
The white disk is the region of instability. In this case the relative equilibria 
are stable not only when there is a dominant mass but also when one of 
the masses is considerably smaller than the other two. This shows that one 
cannot extend Moeckel's conjecture to homogeneous potentials. Moreover as 
a decreases the radius of the disk shrinks and the region of stability becomes 
larger. However if a > the case where mi = 1712 = is always unstable. 



4 Quasihomogeneous Potentials 
4.1 General Case 

Consider the potential U{x). The matrix A = D'WU{x) is of the form 



A 



Anl 



A.. 



where Aij is the 2x2 matrix given by 



Ajk = Ajk + Bjk, if j ^ k 
Akk = ~ ^ ^ Ajk 



(21) 
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where 

A, u = a-^[l-{a + 2)u,uu';,], j^k (22) 

and 

B, , = b^[l-{b + 2)u,,u%]. j^k (23) 

As in the case of the homogeneous potential one can apply Proposition [2] 
taking the vectors v = (1, 0, ■ ■ ■ , 1, 0), w = (0, 1, ■ ■ ■ , 0, 1) that are in the 
kernel of A. Again one finds the eigenvalues fi = ±i 

Now consider v = x (i.e. a relative equilibrium for V, W and U) and 
w = Jx. Equation ffTOl) with the homogeneity of the potentials V and W 
gives: 



and 



M-^D\/V{x)v = -{a + l)M~^VV{x) 

= (a + l)M-^iulMx = (a + 1)ujIv, 



M-^D\/W{x)v = -{b + l)M-^\/W{x) 

= (6 + 1)M^^ujIMx = (6 + 1)ujIv. 



(24) 



(25) 



Consequently 

M-^DVU{x)v = M~^[DVV{x)+DVW{x)]v = [{a+l)ujl+{b+l)uj^]v (26) 

so V is an eigenvector of M'^^A with eigenvalue [{a + l)^^^ + (6 + l)^;!]- On 
the other hand, w = (i?~^(0))'x, as in the case of the homogeneous potential. 
Now 

\/U{R-\9)x) = R-\e)\/U{x) 

because of the rotation invariance of the potential. Differentiating at ^ = 
gives: 

M-^DVU{x)w = M-^JVU{x) = -M-^Juj'^Mx = -A (27) 

where uj"^ = uol + uol- Therefore w is an eigenvector with eigenvalue —uo'^. 
Choosing 

71=[{a + 1)ujI + (6 + 1)ujI]/uj^ = 1 + ^ „ ' (28) 
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and ^ = 



— 1 we find that 



a = 2 ^— ^ 29 

and the roots of Q{z) are positive if [aa;^ + 6a;|]/a;^ > 2. Thus we have proved 
the following 

Theorem 2. Any simultaneous relative equilibrium ofV{x), W{x) and U{x) 
is spectrally unstable if [acuf + 6a;|]/a;^ > 2. 



4.2 Lagrangian Triangle Relative Equilibria 

For n = 3 the equilateral triangle configurations with the center of mass at 
the origin are simultaneous relative equilibria for the potentials V{x), W{x) 
and U{x). They are relative equilibria for any choice of the masses mi, m2 

and 

In this case we assume that the particles are at positions Xj — x'j — c 
where c is the center of mass, and 

x[ = (r, 0) ^2 = (-r/2, rV3/2) 4 = (-r/2, -rV3/2) 

where now we cannot fix the size of the configuration since the potential is 
not homogeneous. Consequently the mutual distances are djk — Tq, where 
ro = V^r, the potential energy is f/(a;) = (mim2+mim3+m2m3)(l/rQ + l/rQ) 
and the moment of inertia is x'^Mx = rQ(mim2+mim3 +77227713) / fi. Moreover 
cuf = ail/ (ro)"'*'^, a;| = {roY^"^ and u"^ = ^1 + u^. As in the homogeneous 
case the x^-^ can be used to compute C = u^'^M^^A, that is of the form 

C 

where 



aV'3(a+2)(m3->7i2) 6 V^(b+2) (m3- m2) 

a(2-a)(m2+m3) b(2-b)(m2+m^) 
7^ 



n v;i(,. ■ 2),R2 , bv^(Ii + 2)m2 

TTpl < 

'o '"o 

a(2-a)m2 b(2-t.)m2 
"^0 "^0 
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Cu 


Ci3 


21 




C23 


31 


C32 


C33 



Cii . 



1 

4^ 



j(2 + 3a)(m2+.T.3) 
loT^ 



av/3(a+2)(m3-tn2) 
oTT^ 



b(2 + 3b)(m2+Tn3) 
3+2 

i.v/3(i.+2)(m3-m2 
3+2 
' n 



C12 



1 



a(2 + 3a)rr 
,a + 2 



iV^(a + 2)n 



b + 2 

fcV^(i»+2)n 

3+2~ 
' n 
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Cl3 : 



1 

4ZJ 



a(2 + 3a)m3 _ b(2 + 3ti)m3 

I^/^3(a + 2)rT^3 _ b b+ 2) ^3 
oT^ b + 2 



iN/3(a + 2)m3 _ bx/3(b + 2),Ti3 
a(2-a)m3 b(2-^b)m3 

ST^ ' bT2 — 



C2I — 



a(2 + 3a)mi _ t(2 + 3b)mi 
a + 2 fa + 2 

' ""0 

iV3(a+2)mi b^{b+2)mi 

— ' b+2 — 



zV3(a + 2)mi bV3(b+2)mi 
a + 2 b + 2 



b(2-bjT] 



C22 — 



a((2 + 3a)Tixi^4m3) 6 ( (2 + 3b)m 1 -47,13 ) 



(iV^(a+2)mi !)^/3(b+2)r 
b+^i 



[iV^(a+2)mi b%/3(b + 2)mi 

ET2 



■0 '0 

(a(a-2)mi+4(l + a)„t3) ( b( b - 2) m 1 +4( 1 + b) ^3 ) 
ET^ ' ET^ 



C23 



1 



4a?Ti3 ^ 4bTTi'^ 
'"6 



4a(l + a)nt3 46(l + b)m3 

S^F^ UT2 



C3I = 



(2 + 3a)mi _ b(2 + 3b)mi 
b + 2 



C32 — 



' 

aV3(a + 2) 

— z^n^ 



aV3(a+2)r 
a + 2 



b-/3{b + 2)rai 
,fa + 2 



bV3{b+2)r 

FTT- 



1 



4,am2 , 4bm2 

a+i ^ 

""0 ""O 



a(2-a)nii b(2-b)mi 

JT^ ' ET2 

""0 ""0 



4a(l + a)ni2 _ 4b(l + b)ni2 
a^F^ b + 2 



■i((2 + 3a)mi -4m2) , b((2 + 3b)77ii -4^3) a.V3{o. + 2)m.-^ bV3{b+2)mi 

oT^ ^ b + 2 b+2 

^0 ^0 ^0 ^0 

a-/3{cL+2)m.i , by3(b+2)Tni a ( (a - 2)n. 1 +4( 1 + a) ^2 ) , b( (b- 2)n. 1 +4( 1 + b) ^13 ) 

ETl ' ET2 ETl ' ET2 

^0 '0 ^0 ^0 



Four eigenvalues of C were found in subsection 4.1: 0, 0, —1 and 1 + 2:i£i^!£2._ 
Or, using the values of cui, ^2 and uj found in the case of three bodies 0, 0, —1 
and 

Theorem [2] applied to the triangle solutions shows that they are spectrally 
unstable for any value of the masses mi, 1712 and when 

(a^ - 2a)r^+2 + (6^ - 2b)r^ ^ ^ ^^^^ 



ar, 



Let us denote the remaining eigenvalues of C with t] and ^. We will only 
need to know the value of 77 + ^ and fjC,. Taking the trace of C gives 

„2 6+2 , Z,2 a+2 /,2 fe+2 , /,2 a+2 

Tr(C) = °l,+ °. +? + ?=2-°'° +"'» 



ar. 
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or 

In the following we will use 6 to denote the right hand side of the previous 
equation. The first parameter in Proposition [2] is 



To find /5 = (r/ + l)(^ + l) = 1 + 5 + rj^ one needs to compute r/^. In order to 
do that one can observe that the sum of all possible pairs of eigenvalues of 
C can be found in two ways. Direct computations give —1 — 5 + 5^ + V^- On 
the other hand one can find this sum as T = (l/2)[Tr(C)^ — Tr((C)^)] where 

T = [{ml + ml + ml){rl''+\4b^ - 86^ - 8b^) + r^'-^^Aa" - 8a^ - 8a^)) 
+ {ml + ml + ml){r^+''+\8a%^ - 8ab^ - 8a% - IQab))) 
+ {mim2 + m2m^ + mim^) 

X {rl'^+^llb^ - 46^ - 46^) + rl'^+^lla^ - Aa"^ - Aa^)) 
+ {mim2 + mamg + mim3)(ro+^+'^(22a^6^ - 4a6^ - 4a^6 - 8ab))] 
''o ) ) 

(31) 



X (4/i2(ar^+2^aro"+2^2, 



which leads to /5 = 2(1 + 5) — 5^ + T, or explicitly, after some simplifications 
to 



_ 3(mim2 + mim3 + m2m3) / b{b + 2)rg+^ + a(a + 2)ro+^ ^ 



Since /? > the equilateral triangle relative equihbrium is spectrally stable 
if and only if a > and 

m^m^ + m^m-i + m2m-i ^ 1 (bib - 2)ro"'' + a(a - 2)\ ^ 



(mi + ms + ma)^ - 3 V &(6 + 2)r;j-'' + a(a + 2 

This generalizes Routh stability condition to the case of quasihomogeneous 
potentials. 

We now want to show that, in this case, the stability depends, for certain 
values of the masses, on the size of the configuration, i.e. of the equilateral 
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triangle. In order to do that let /(mi, 7712,^3) be the left hand side and 
g{a,b,rQ) be the right hand side of equation (15^ and consider the hmits of 
g{a, b, To), as ro — > and tq — * 00. The first limit is 



With the preparations above we are well on our way to proving the following 

Theorem 3. Consider the Lagrange triangle solution for the quasihomoge- 
neous potential U. 

(a) If Q <b < a <2 it is 

(i) Unstable when f > l/3((6 - 2)/{b + 2))^ 
(a) Spectrally stable when f < l/3((a — 2)/(a + 2))^ 
(ill) Unstable when l/3((a-2)/(a+2))2 < / < l/3((6-2)/(6+2))2 and 
r < rl, spectrally stable when l/3((a-2)/(a + 2))^ < / < l/3((&- 
2)/(6 + 2))^ and r > Vq for some Vq e (where = rl{f) ). 

(b) If0<b<2<a it IS 

(i) Unstable when f > (1/3) [(6 - 2)/{b + 2)]^ 

(11) Unstable when < / < (1/3) [(6 - 2)/{b + 2)]^ and tq < zl for 
some zl e (where zl = zKf) ). 

(ill) Spectrally stable when < / < {l/?>)\{b-2)/{b + 2)f andro > zl- 

(c) if 2 < b < a it is unstable for any value of the masses. 

Proof. First we prove part (a) of the theorem. When 0<6<a<2, «>0 
then \imrQ-,o g{a, b, vq) < limj-g^oo 5'(o, b, tq), since 1/3 ((a — 2) /(a + 2))^ is a 
monotonically decreasing function for < a < 2. Therefore from inequality 
fl33|) follows that when the masses satisfy 




while the second is 




1 /a-2 
3 \a + 2 



) 



< 



mim2 + mim^ + 7712777.3 



(7?7i + 7772 + 



< - 




) 
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Figure 2: (a) Regions of stability in the triangle of normalized masses for the 
potential U with 0<6<a<14 — 8-\/3. The equilateral triangle is unstable 
for the mass triplets in Di and spectrally stable for the ones in D^. When 
the mass triplet is in D2 the stability depends on the size of the triangle: it 
is unstable for r < Tq and stable for r > Tq. (b) Regions of stability in the 
triangle of normalized masses for the potential U when < 6 < 14 — 8\/3, 
a > 2. The equilateral triangle is unstable for the mass triplets in Ei. When 
the mass triplet is in E2 the stability depends on the size of the triangle: it 
is unstable for r < and stable for r > 



the relative equilibrium is spectrally unstable as ro 
when To ^ 00. Moreover 

dg _2b{a- b)r^-''[{b'^ - Ihy-^' + [0? - 
dro 3 



but spectrally stable 



2a)] 



b{b + 2)r^-' + a{a + 2) 



X 



(6 -2) + (6 + 2) 



(6^ - 2b)r 



a—b 



2a) 



bib + 2)rl-' + a(a + 2) 



(34) 



is always negative if a < 2, 6 < 2 and a < b. Therefore g{rQ) is a monotoni- 
cally decreasing function of Vq. Consequently there exists an unique Tq such 
that the triangle is unstable for r < Tq and spectrally stable for r > r^. In 
other words the relative equilibrium is unstable when the equilateral triangle 
is small but it is spectrally stable if the triangle is large. 

On the other hand, since g{ro) is a monotonically decreasing function of 
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ro then l/3((a - 2)/(a + 2)f <g< l/3((6 - 2)/(6 + 2)f. Thus the triangle 
is spectrally stable when / < l/3((a — 2)/(a + 2))^ and it is unstable when 



/ > l/3((& — 2)/(6 + 2))^. Figure 2(a) depicts the regions of stability and 
instability for < 6 < a < 14 - 8^/37 

We now prove part (b) of the theorem. Observe that a > in this case 
gives ro > z* where 

— 2a' 
' 62 _ 26 

Thus the triangle is unstable for any value of the masses when < z*. 

On the other hand g{a,b,z*) = 0. Moreover from equation (!34l) one can 
see that dg/dro is negative for ro < z* and positive for tq > z*. This shows 
that (yf(a,6, ro) is a monotonically decreasing function of ro in the interval 
(0, z*) and a monotonically increasing one in {z*, oo). Therefore if we fix the 
values of the masses so that < / < (l/3)[(6 — 2)/(6 + 2)]^ the equation 
/ = g{a, b, ro) has a unique solution in the interval {z*, oo). Consequently 
/ > g{a, b, To) for ro < z^ and the triangle is unstable, while / < g{a, b, vq) 



for ro > zl and the triangle is spectrally stable. Figure 2(b) depicts the 
regions of stability and instability for < 6 < 14 — 8y/3, a > 2. 

The proof of part (c) follows immediately from inequality ( l30l) . □ 
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